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Abstract 

The paper is devoted to a study of the peristaltic motion of blood in the micro- 
circulatory system. The vessel is considered to be of varying cross-section. The progressive 
peristaltic waves are taken to be of sinusoidal nature. Blood is considered to be a Herschel- 
Bulkley fluid. Of particular concern here is to investigate the effects of amplitude ratio, 
mean pressure gradient, yield stress and the power law index on the velocity distribution, 
streamline pattern and wall shear stress. On the basis of the derived analytical expression, 
extensive numerical calculations have been made. The study reveals that velocity of blood 
and wall shear stress are appreciably affected due to the non-uniform geometry of blood 
vessels. They are also highly sensitive to the magnitude of the amplitude ratio and the 
value of the fluid index. 
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1 Introduction 

Peristaltic transport is well known to physiologists as a natural mechanism of pumping materials 
in the case of most physiological fluids. Apart from physiological fluids, some other fluids 
also exhibit peristaltic behaviour. Peristalsis usually occurs when the flow is induced by a 
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progressive wave of area contraction/expansion along the length of the boundary of a fluid-filled 
distensible tube. Usefulness of studies on peristaltic flow has been discussed in detail in our 
earlier publications (Misra et al. [DE1E1I3E1EJC3E1E], Maiti and Misra [TOlITT]) and in some 
other references [121 02] • 



Nomenclature 






Half-width of the channel at the inlet 


b 


Wave amplitude 


H 


Vertical displacement of the wall 


n 


Flow index number 


m 


Reciprocal of n 


P 


Fluid pressure 


Q 


Flux at axial location 


t 


Time 


X,Y 


Rectangular Cartesian co-ordinates 


U,V 


Velocity components in X,Y directions respectively 


5 


Wave number 


Ap 


Pressure difference between the channel ends 


A 


Wave length of the travelling wave motion of the wall 


A 


A parametric constant 




Blood viscosity 


V 


Kinematic viscosity of blood 





Amplitude ratio 


P 


Density of blood 




Yield stress of blood 


Th 


Wall shear stress 



The phenomenon of peristalsis plays an important role in the functioning of heart-lung machine, 
blood pump machine and dialysis machine. Fung and Yih [H] presented a theoretical analysis of 
peristaltic transport primarily with inertia-free Newtonian flows driven by sinusoidal transverse 
waves of small amplitude. Investigation of peristaltic motion in connection with functions of 
different physiological systems such as the ureter, the gastro-intestinal tract, the small blood 
vessels and other glandular ducts was first made by Shapiro et al. [15]. They presented a closed 
form solution for an infinite train of waves for small Reynolds number flow. Their study was, 
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however, restricted to cases where the wave length is long and the wave amplitude is arbitrary. 
Conditions for the the presence of physiologically significant phenomena of trapping and reflux 
were also suggested by them. Jaffrin and Shapiro [13j as well as Srivastava and Srivastava 
[16] summarized the early literatures on peristaltic transport. Some of the recent studies on 
peristaltic transport were discussed by Usha and Rao [17], Mishra and Rao [18], Yaniv et al. 
[19], Jimenez- Lozano et al [20], Wang et al. [21] and Hayat et al. [22]. 

Some important theoretical analyses on different aspects of blood flow were carried out in 
a systematic manner by Misra et al. [23j EH [25] . Attempt to consider the complex rheology of 
various physiological fluids was made in several studies [261 1211 [29J, [30j [311 [32j [33j [M] . The 
non- Newtonian behaviour of blood mainly owes to the presence of erythrocytes in whole blood. 
In the case of blood, such behaviour starts becoming prominent when the hematocrit rises above 
20%. This particular behaviour plays a dominating role when the hematocrit level lies between 
40% and 70% [351 ESI [37J . Some other relevant theoretical studies on non-Newtonian fluid flows 
were carried out by Masud and Kwack [38], Kwack and Masud [39] as well as by Anand and 
Rajagopal [3D] . 

It is known that the flow behaviour of blood in small vessels (diameter<0.02 cm) and at low 
shear rate (< 20sec _1 ) can be represented by a power law fluid [411 H2]. Merill et al. [43] pointed 
out that Casson model holds satisfactory for blood flowing in tubes of 130-1000/im. Moreover, 
Blair and Spanner [44] reported that blood obeys Casson model for moderate shear rate flows. 
However, they pointed out that for cow's blood, Herschel-Bulkley model is more appropriate 
than Casson model. 

It is known that physiological organs are by and large non-uniform ducts [45 l 146 1 147] . Several 
authors [lEJ HHJ HH] made some initial attempts to perform theoretical studies pertaining to 
peristaltic transport of physiological fluids in vessels of non-uniform cross section. These analyses 
were mostly restricted to the assumption of either a Newtonian fluid or a non-Newtonian fluid 
of Casson/power-law type. Moreover, in these reports the different flow characteristics have not 
been adequately discussed. The strong merit of Herschel-Bulkley model is that fluids represented 
by this model describe very well material flows with a non-linear constitutive relation depicting 
the behaviour of shear-thinning/shear-thickening fluids that are of much importance in the field 
of biomedical engineering [50] . It is also to be noted that for formulating a non-Newtonian 
model of blood, Herschel-Bulkley fluid model is more general than most other non-Newtonian 
fluid models. It is worthwhile to mention that results for a fluid represented by Bingham plastic 
model, power law model and Newtonian fluid model can be derived from those of the Herschel- 
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Bulkley fluid model. Also this model has been found to yield more accurate results than many 
other non-Newtonian models. 

In view of the above, we have taken up here a study on the peristaltic transport of blood in a 
non-uniform channel, by treating blood as a Herschel-Bulkley fluid. It is worthwhile to mention 
that flow through axisymmetric tubes is qualitatively similar to the case of flow in channels. 
On the basis of their theoretical/experimental study on peristaltic pumping at low Reynolds 
number, Shapiro et al. [15] also confirmed that flow behaviour in the case of an axisymmetric 
tube is identical to that in the case of channel flow. The formulation and analysis presented 
in the sequel are particularly suitable for investigating the peristaltic motion in vessels of small 
dimensions, e.g. arterioles and venules. Since for flow of blood through smaller vessels in the 
micro-circulatory system, the Reynolds number is low and since the ratio between half-width of 
the channel under consideration and the wave length is considered small, it has been possible 
to perform the theoretical analysis for the present problem in a convenient manner, by using 
the lubrication theory [15J. On the basis of the derived analytical expressions, computational 
work has been executed keeping a specific situation of micro-circulation in view, with the main 
purpose of examining the distribution of velocity of blood and wall shear stress as well as the 
streamline pattern, trajectories of individual fluid particles and pumping performance. The plots 
for the computed results give us a clear idea of qualitative variation of various fluid dynamical 
parameters. The results indicate plug flow in the central region, where most of the erythrocytes 
are accumulated and non-plug flow in the peripheral region. The results of the present study 
are in good agreement with those reported earlier by previous investigators. The study bears 
the potential to explore some important phenomena that are useful for having a better insight 
into the flow in the micro-circulatory system. 

The study has an important bearing on the clinical procedure of extra-corporeal circulation 
of blood by using the heart-lung machine, where there is a chance of damage of erythrocytes 
owing to significant variation of the wall shear stress. The results are also likely to find important 
application in roller pumps and arthro-pumps by which fluids can be transported in living organs 
in pathological states. 

2 Formulation and Analysis 

Let us consider the peristaltic motion of blood, by treating it as an incompressible viscous non- 
Newtonian fluid. The non-Newtonian behaviour is considered to be of Herschel-Bulkley type. 
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Y=0 




Figure 1: A physical sketch of the problem for a tapered channel 

We shall study two-dimensional channel flow, width of the channel being non-uniform. 

We take (X,Y) as Cartesian coordinates of the location of a fluid particle, X being measured 
in the direction of wave propagation and Y in the normal direction. Let Y = H and Y = —H 
be respectively the upper and lower boundaries of the channel (cf. Fig. [1]). The medium is 
considered to be induced by a progressive sinusoidal wave train propagating with a constant 
speed c along the channel wall, such that H = d(X) + a sin(^(X — ct)), with d{X)=d + AX. 
d(X) represents the half-width of the channel at any axial distance X from the inlet, do being 
the half- width at the inlet, A(<1) a constant whose magnitude depends on the length of the 
channel as well as the inlet and outlet dimensions, a the wave amplitude, t the time and A the 
wave length. 

Under the assumptions stated above, the basic equations that govern the fluid motion are 
the field equations 

V • V = 

, dV ^ 

and p— — = v-a + pt, 
dt 

where V is the velocity, f the body force per unit mass and p the density of the fluid, while 4 
denotes the material time derivative, a is the Cauchy stress defined by o = —PI + T, where 
T = 2pEij + r]IS and S = V • V, being the symmetric part of the velocity gradient L, 
Eij = \[L + L T l L = VV. 

—PI denotes the indeterminate part of the stress due to the constraint of incompressibility; 
p and f] are viscosity parameters. The Herschel-Bulkley model gives the combined effect of 
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Bingham plastic and power-law behavior of a fluid. At low strain rate ( 7 < ^ ), the material 
acts like a viscous fluid with constant viscosity /i . But when the strain rate increases and the 
yield stress threshold, r , is reached, the fluid behavior is described by a power law 



7 

where a and n denote respectively the consistency factors and the power law index, n < 1 
and n > 1 correspond to a shear thinning fluid and a shear thickening fluid respectively. If 
the channel length is an integral multiple of the wavelength, the pressure difference across the 
ends of the channel is a constant. The pressure p remains constant across any axial station 
of the channel, when the wavelength is large and curvature effects are negligibly small. Since 
in the study, the geometry of wall surface is non-uniform, the flow is inherently unsteady in 
the laboratory frame as well as in the wave frame of reference. Disregarding the body forces 
(i.e. taking f=0), and using the Herschel-Bulkley equations, the governing equations of the 
incompressible fluid motion in the micro- vessel in the fixed frame of reference may be put in the 
form 

fdU MTT dU dU\ dP dr xx dr XY 
(dV dV dV\ dP dr YX dr YY 

fdUi dUj\ 

with nj = 2fiE tj = fi ^— + 5x J ' ( 3 ) 

{ no : for n < n , 

A* — i (4) 
[ aW 1 - 1 + Ton" 1 : for U > U 

IT = yf2E~E~ 3 (5) 

The limiting viscosity fx is considered such that 

fi = aWJ- 1 + tqUq 1 (6) 

In the analysis that follows, we shall make use of the following non-dimensional variables: 

<y x v> Y rr U - V r d - dl +l P - ct , H , a 

X = T' Y = T^ U= -> V = ~^ 5 = T' P = J 4r' t = T' h = T> ^=T> 
A do c co A /ic n A A do d 

TO = 7 Hf, T Y X = 7 wi, V = — (7J 
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The equation governing the flow of the fluid can now be written in the form (dropping the bars 
over the symbols) 



r) ,'dU TT dU „dU 
Re6i lH +U dX +V dY 



dp hM , oM) | gft(S+^S)) 



dX 



dX 



dY 



r) ,'dV TT dV X 8V 

Re5i l* +U dx +v W 



~dY +6 dX + 6 ^Y~ (9) 



where <3> 



\ 



fdU\ 2 fdV x2 



25U ^ 



- 



\dX J \dY 



fdU „dV\ 



- 



\dY dX 



n-l 



+T 



\ 



25 2 



dU 
dX 



dVY 

WJ 



dU 
dY 



Using the long wavelength approximation (5 <C 1) and the lubrication approach [15 
governing equations and the boundary conditions describing the flow in the fixed 
reference may be rewritten in terms of the dimensionless variables as 

dP _ dr YX 
dX ~ dY ' 



(10) 

HE], the 
frame of 

(11) 



^ = 
dY 



where tyx 



To 



dU 



dY 



dU 

s 9 n ( ^7 ) (c/. M 



(12) 
(13) 



* = 0, U Y = ^yy = 0, t yx = at Y = 0; U = ^ Y = at Y = h 



(14) 



Thus the pressure gradient |j is independent of Y. The solution of equation (iTTj) satisfying (fl4|) 
is found in the form 



[{P x h - r )^ +1 - (PxF - ro)"^ 1 ] : i/ Y > 



-PiY - T 



ini+11 



: if Y<0, 



(15) 



where Pi — — |j, ni = ~. If the plug flow region be given by Y = Y , 



U Y = at Y = Y 
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Then Y r 



Pi 



h denned in (0) stands for the non-dimensional vertical displacement. If tyx = Th at Y=h, we 
find h = Th/P\. 



Now — = — = t (say), < r < 1 
h r h 



The plug velocity is then given by 

{P x h - r )™ 1+1 



U„ 



Using the boundary conditions 



(ni + l)Pi 



^ p = at Y = 
and * = * p at Y = Y , 
and integrating ( JT51) and ( ITTj) . the stream function \l/ is found to be given by 



(ni+l) 



(ni+1) 



n! +2 
+2 

!+i , (^y-T )"i+ 2 



y (/i - y ) 



if Y <Y<h 

if - h < Y < -Y , 



(16) 



(17) 



and ty. 



P^/i- Y ) ni+1 Y 



if 



Y n < Y < Y n 



ni + 1 

The instantaneous rate of volume flow through each section, Q(x,t), is given by 

Q{X,t)= [°U Y dY+ ! U dY 
Jo Jy 

p™ i (/i-y ) ni+1 (n 1 / i + /i + y ) 



(n 1 + l)(n 1 + 2) 



m = - 

n 



(19) 



(20) 



Now 



dP_ _ f Q(X,t)(n 1 + l)(n 1 + 2) 

ax ~ [(A - y )™i +1 (ni/;, + /i + y )_ 

Q(X,t)(n 1 + l)(n 1 + 2) 
_/i ni+2 (l - r) n i +1 (ni + 1 + t) 

The average pressure rise per wave length is calculated as 



AP 



(21) 



(22) 
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In the fixed frame of reference, the expression for the non-dimensional transverse velocity V is 
found to be given by 

ni+2' 



V(X,Y,t) 



m (m + i)(m + 2) \ \ UH ly _ h 2 (Y_ _ ^ 



Q n (X,t) h(l-T)( ni + l + T) J v ' ni + 2\h 



+(», + 1)(1 - r)?y] x ggM^Wizf) (23 ) 



It may be noted that if we put n = 1, r = 0, A = in equations ( 1T5|) . (JTHJ), ( 120]) and ( l2~il . 
the expressions reduce to those reported earlier in [15]. Our results also tally with those of [52J, 
when the eccentricity of the elliptical motion of cilia tips is set equal to zero in their analysis 
for a Newtonian fluid flowing through a uniform channel. Moreover, when n — 1 and Y = 0, 
the expression for the pressure gradient given by ([2T1) reduces to that obtained by Gupta and 
Seshadri [18] for the peristaltic motion of a Newtonian fluid having constant viscosity. It may 
be noted that since the right hand side of equation (122]) cannot integrated in closed form, for 
non- uniform/uniform geometry, for further investigation of our problem, we had to resort to the 
use of appropriate softwares, as mentioned in the next section. 



3 Computational Results and Discussion 

The instantaneous rate of volume flow Q(X,t) has been assumed to be periodic in (X-t) ([16] 
HHlllS]), so that Q appearing equation (122]) can be expressed as 

Q n (X, t) = Q n + (j) sin 2tt (X - t) , (24) 

Q being the time-averaged flow flux. Due to complexity of the problem, it has not been possible 
to find the expression for the average pressure rise, AP given by ( 12T]) . It has been computed 
numerically by using the software Mathematica. 

In this section, on the basis of the present study, we shall obtain theoretical estimates of 
different physical quantities that are of relevance to the physiological problem of blood flow 
in micro-circulatory system. For this purpose, we have used the following data valid in the 
physiological range ([I2J EH ESI El]): d = 10 to 60/im, = 0.1 to 0.9, f = 0.01 to 0.02, 
AP = —300 to 50, r = 0.0 to 0.2, Q = to 2, n = | to 2. Unlike other studies on non-uniform 
geometry, the value of A has been so chosen that for converging tubes (e.g. arterioles), the width 
of the outlet of one wave length is 25% less than that of the inlet, while in the case of diverging 
tubes (e.g. venules), the width of the outlet of one wave length is 25% more than that of inlet. 
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<P 


Our results 


Kesuits m [zoj 


Wave Crest 


O.z 


n O A A 


O.ZO 


0.3 


0.338 


0.344 


0.4 


0.418 


0.427 


Wave Trough 


0.2 


-0.366 


-0.375 


0.3 


-0.627 


-0.635 


0.4 


-0.976 


-0.984 



Table 1: Comparison of the present study with those of Takabatake and Ayukawa [25] for small 
values of Re and S (n — 1, r = 0, A = 0, Q — 0) 

It is important to mention that the results presented in the sequel for shear thinning and 
shear thickening fluids are quite relevant for the study of blood rheology. Normal blood usually 
behaves like a shear thinning fluid for which with the increase in shear rate, the viscosity 
decreases. It has been mentioned in [5U [551 ES] that in the case of hardened red blood cell 
suspension, the fluid behaviour is that of a shear thickening fluid for which with the fluid 
viscosity is enhanced due to an elevation of the shear rate. 

3.1 Velocity Distribution 

Plots in Figs. [2J|6]give the distribution of axial velocity in the cases of free pumping, pumping and 
co-pumping for different values of the amplitude ratio 0, flow index number n, r, A. Fig. [2]^a) 
shows that the results computed on the basis of our study for the particular case of Newtonian 
fluid tally well with the results reported by Shapiro et al. [15] when the amplitude ratio = 0.5. 
Variation of axial velocity in the vertical direction at the wave crest and the wave trough 
are exhibited in Figs. [2]^b,c) for different values of 0. In order to compare our results with 
those reported by Takabatake and Ayukawa [27] , we have reproduced their results alongside the 
results computed on the basis of the present study in Table 1. From the tabulated values, one 
can observe that axial velocity along the central line at the wave crest and the wave trough 
match well with [27] when the wave number is small. Since the velocity profiles and height of 
the channel change with time, investigation has been made on the basis of the present study 
for the distribution of velocity at an interval of T/4. Figs. |3] depict the aerial view of a few 
typical axial velocity distributions for a Newtonian fluid flowing over a uniform channel. The 
corresponding velocity contours for the peristaltic motion corresponding to = 0.1 are shown in 
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0=0.2 — 




o 0.3 -- 




0=0.3 







0.5 




Y 



0.5 



1.5 



(b) 



(c) 

Figure 2: Variation of axial velocity in the vertical direction (a) at x=0.5, t=0.25 (A = 0, = 
0.5, r = 0, n — 1, AP = 0) (b) for different values of <fi at wave crest when A = 0, n — 1, 
r = 0, Q = (c) for different values of at wave trough when A = 0, n = 1, r = 0, Q = 
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(a) for t = 0.0 



(b) for t = 0.25 




(c) fort = 0.5 (d) fort = 0.75 

Figure 3: Aerial view of the velocity distribution at different instants of time (n = 1, A = 
0, AP = 0, r = 0, = 0.5) 
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Figs. HJ The velocity contours presented in Fig. HJ^b) resemble with those in Fig. H)^e) reported 
by Selvarajan et al. [57], although in our plots, there is a slight increase in area of contour lines 
at the region of wave expansion. This is possibly because the curved boundary was not taken 
into consideration in [57] for calculating the velocity contours. 

Fig [5] gives the velocity distribution in the plane of the channel. This figure reveals that 
at any instant of time, although there exists a retrograde flow region, the forward flow region 
is predominant, the time averaged flow rate being positive. We have also found that if Q = 0, 
the retrograde flow region occupies exactly half of the one complete wave length, the remaining 
half being occupied by the forward mean flow region. For t = T/4, our study reveals that there 
exists two stagnation points on the axis, one being located near the trailing end and the other 
near the leading end. Both these observations agree with those of Takabatake and Ayukawa 
[27], who performed a similar study numerically for a Newtonian fluid. Fig. E(a) shows that 
when n = 1 (that is, in the case of a Newtonian fluid) and AP = 0, in both the forward flow 
region and the backward flow region velocity increases as the value of (f) is raised. One can also 
have an idea of the extent by which the value of <fi affects the axial velocity for shear thinning 
(cf. Fig. [6(b) ) and shear thickening fluids (Fig. Efc)). It is apparent from these figures that in 
either region, as r increases, the magnitude of velocity decreases for both types of fluids. 

The extent to which the velocity distribution is influenced by the value of the rheological 
fluid index 'n' can be observed from Figs. E(d-i) for uniform/non-uniform channels. We observe 
that the the parabolic nature of the velocity profiles is disturbed due to non-Newtonian effect 
and that as the value of 'n' increases, the magnitude of the velocity increases. It is worthwhile 
to mention that for a converging channel, the magnitude of the velocity is greater than that 
in the case of a uniform channel; however, for a diverging channel, we have an altogether 
different observation. Figs. E£j-1) depict the influence of pressure on velocity distribution for 
shear thinning (n < l)/shear thickening (n > 1) fluids. While the plots given in Figs. GJEO) 
correspond to the case of AP = (free pumping), Figs. E^j-1) have been plotted for the co- 
pumping case AP < 0. For these plots, we consider the case when AP = —1. In this case, for 
a shear thinning fluid, flow reversal is totally absent in the case of a uniform/diverging channel. 
For a converging channel, however, although there is reduction in the region of flow reversal, 
it does not vanish altogether, irrespective of whether the fluid is of shear thinning or of shear 
thickening type. 
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(b) for t = 0.25 




(c) for t = 0.5 
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(e) Contours reproduced from Selvarajan et al. [55] 



Figure 4: Axial velocity contour at different time when = 0.1, n = 1, A = 0, r = 0, Q = 
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Figure 5: Velocity distribution at different instants of time 
3.2 Pumping Performance 

Since peristaltic transport of a fluid is associated with the concept of mechanical pumping, it 
is worthwhile to examine the pumping performance under the purview of the present study. 
The average pumping flow rate corresponding to a given pressure against which the pumping 
action takes place gives a measure of peristaltic pumping performance. By using lubrication 
theory, it has been observed by Shapiro et al. [15] that in the infinite tube model, the flow 
rate averaged over one wave varies linearly with pressure difference. However, the existence 
of a peripheral layer (Newtonian fluid), makes the relationship non-linear [58]. The present 
study being concerned with the peristaltic transport of a non-Newtonian fluid, the relationship 
between the pressure difference and the mean flow rate is found to be non-linear (cf. Figs |7(b- 
h)), although the study is based upon the consideration that the channel length is an integral 
multiple of the wave length. This has been the observation both for converging and diverging 
channels for all values of the fluid index n 1). In the Newtonian case, however, the relation is 
linear (cf. [7(a)) and this is in conformity to the observation reported by previous investigators 
[23 [59]. 

The plots in Figs. UJ show that the mean flow rate, Q increases as AP decreases. Figs. UJb) 
and[JJ(c) reveal that the pumping region increases with the increase in the value of the amplitude 
ratio <p for both shear thinning and shear thickening fluids. Figs. UJd-f) depict the influence 
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Figure 6: Distribution of velocity in different situations 
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Figure 7: Pressure rise versus flow rate 
18 



of the rheological parameter 'n' on the pumping performance in uniform/ diverging/ converging 
channels. It is important to note that for all types of channels, the pumping region (AP > 0) 
significantly increases with the increase in the value of 'n'. In the co-pumping region (AP < 0) 
the pressure rise decreases when Q exceeds a certain value. From Figs. [7(g) and[7Jh), we find 
that Q is not significantly affected by r for free pumping case. Moreover, when Q exceeds a 
certain critical limit (equal to the value of Q in the free pumping case), pressure rise decreases 
with the increase in r. We further find that for both shear thinning and shear thickening fluids, 
pumping region increases with r increasing; however, the effect of r is greater in the case of 
shear thickening fluids. 

3.3 Wall Shear Stress 

It is known that when the shear stress generated on the wall of a blood vessel exceeds a certain 
limit, there is a possibility that it can cause damage to the constituents of blood. Moreover, 
the magnitude of wall shear stress plays an important role in the molecular convective process 
at high Prandtl number/Schmidt number |60j. Because of this reason, investigation of shearing 
stress deserves special attention in the hemodynamical flow of blood in arteries. Figs. [8] depict 
the wall shear stress distribution under varied conditions. Fig. [E(a) gives the distribution of 
wall shear stress at four specially chosen instants of time during one complete wave period. This 
figure shows that at each of these time instants, there exist two peaks in the wall shear stress 
distribution, with a gradual ramp in between; however, negative peak of wall shear stress, r min 
is not as large as the maximum wall shear stress, r max . It may be observed that the transition 
from T m i n to T max takes place in the zone between the maximal height and the minimum height 
of the channel (cf. Fig. [8(a)). At the point of maximum occlusion, the wall shear stress as well 
as the pressure is maximum. The pressure gradient to the left of this point is positive and hence 
the local instantaneous flow will take place towards the left of r max . This may be responsible 
for a number of consequences. For example, if the rate of shear at the crest is quite high, a 
dissolving wavy wall will have a tendency to level out. Also, some chemical reaction between 
the wall material and the constituents of blood is likely to set in. Owing to the deposition of 
the products of the chemical reaction , wall amplitude increases at a rapid rate. This may lead 
to clogging of the blood vessel. 

In other regions, for wall shear stress distribution curves, the peaks on both sides of r max 
are small. Thus the local instantaneous flow will occur in the direction of the peristaltic wave 
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Figure 8: Wall shear stress distribution for different cases 
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Figure 9: Streamline patterns for peristaltic flow of a Newtonian fluid at different instants of 
time (Q = 0, A = 0, r = 0, <p = 0.5). 
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igure 10: Streamline patterns in the case of peristaltic flow for different values of <fi when 
= 1, t = 0.25 ,r = 0, A = 0, Q = 



23 



propagation. Of course, when averaged over one wave period, net effect may be looked upon as 
peristaltic transport in the direction of wave propagation. It may be observed from Figs. [Sfb-c) 
that in the contracting region where occlusion takes place, there is a significant increase in the 
wall shear stress due to an increase in the value of <j). This applies for both shear thinning and 
shear thickening cases. However, in the expanding region, with the increase in 0, r m j ra increases 
in magnitude for a shear thinning fluid; for a shear thickening fluid, the effect is relatively little. 
Fig. Myd) shows that with increase in r, but for r m i n the effect is minimal. The 

effect of the rheological fluid index 'n' on the distribution of wall shear stress has been shown 
in Figs. IHe-g) for uniform/non- uniform channels. It may be noted that in all types of channels 
studied here, values of T max enhance with an increase in the value of the fluid index n. Moreover, 
one may observe that the shear stress difference between the outlet and the inlet in the case 
of converging channel is exceedingly large in contrast to the case of a diverging channel. Figs. 
Etff-j) indicate that as the time averaged flow rate increases, the wall shear stress decreases. 
This is so for all the cases examined here. It is also to be noted that the observation regarding 
the shear stress difference between the outlet and the inlet for Q > is similar to that in case 
Q = discussed earlier. 

3.4 Streamlines 

Figs. E31IH gives an insight into the changes in the pattern of streamlines that occur due to 
changes in the values of various parameters that govern the flow of blood under the purview of 
the present study. Streamlines for different flow governing parameters are depicted in Figs. I9HT41 
These figures indicate that in the portion of the channel where it is dilating, the flow is pulled 
by the wall, where as in the contracting portion, the flow is pushed away from the wall. It may 
be noted that since the flow behaviour is unsteady in the fixed frame of reference, at different 
times, streamlines are of different nature. Typical nature of streamlines for the problem under 
the present consideration at different instants of time is shown in Figs. |9j 

In the wave frame of reference, the formation of an internally circulating bolus of fluid 
that moves along with the same speed as that of the wave is a very interesting phenomenon 
from the view point of fluid dynamics. This physical phenomenon is usually referred to as 
'trapping'. However, in the fixed frame of reference it does not appear. It may be mentioned 
that investigation of the streamline patterns is quite important, particularly for some type of 
problems, because of the fact that the difference between the values of the stream function at any 
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Figure 11: Streamline patterns in the case of peristaltic flow of rheo logical fluid for different 
values of r when Q = 0, t = 0.25, A = 0, = 0.5 
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Figure 12: Streamline patterns in the case of peristaltic flow for different values of n when Q=0, 
t=0.25, r = 0.1, A = 0, = 0.5 
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two points can be used to calculate the volumetric flow rate/flow flux through a line connecting 
the two points. 



3.5 Trajectory of Particles and Reflux 

In the micro-circulatory system, the nature of trajectories of fluid particles play an important 
role in the functioning of arterioles. The reason is that bio-chemical reactions that take place 
between blood and the vessel constituents are effected by convective transport of the fluid 
particles. Owing to the said bio-chemical reactions there may be a fast increase in the wall 
amplitude. This may lead to clogging of blood. From the plots of the velocity distribution (cf. 
Fig. [5]) it is clear that at a given cross-section, flow of blood takes place alternately in/opposite 
to the direction of wave propagation in different phases. The functioning of arterioles may 
then be affected by the confluence of the convective transport in the longitudinal direction. 
By investigating the pathlines of massless particles moving in the direction opposite to that of 
the peristaltic wave propagation in the Lagrangian frame of reference, it is possible to have an 
insight into the reflux phenomenon. 

By resorting to appropriate numerical methods for solving the simultaneous differential 
equations 

£ = I7, ^ = V (25) 
dt dt v ; 

successively starting from the initial location of the particles, it has been possible to deter- 
mine the trajectories of the particles. In the above equations, (X,Y) are the non-dimensional 
coordinates of a particle at time t. 

One of the novel features of this study is the examination of the particle trajectories for a 
non-Newtonian fluid. To the best of our knowledge and observation, no previous investigator 
has ever studied the trajectories for a non-Newtonian fluid. It is very important to observe 
that the convergence of the pathlines is very sensitive even to a moderate deviation from n=l. 
Integration of the above-written differential equations has been carried out by applying the 
Runge-Kutta 4 method. Figs. [ToT a-b) present the trajectories for the Newtonian case, while 
Figs. [T5T c-d) show the trajectories for the non-Newtonian case. The trajectories displayed in 
Fig. [ToTa) are exactly similar to those reported in [15J. It is worthwhile to observe that the 
period of the particle is different from the wave period and that when Q = 0, particles in the 
vicinity of the axis of the channel in most cases undergo a net positive displacement, whereas 
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Figure 13: Streamline patterns for the peristaltic flow of a shear thinning fluid (n=2/3) for 
different values of Q when t=0.25, r = 0.1, A = 0, = 0.5 
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Figure 14: Streamline patterns for the peristaltic flow of a shear thickening fluid (n=4/3) for 
different values of Q when t=0.25, r = 0.1, A = 0, = 0.5 
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displacement of particles in the neighbourhood of the boundary is negative. These observations 
tally with those reported in [151 [27]. However, at the point (0.75,0.2) where the velocity profile 
of the particle at x=0.75 and time t=0.0 is minimum, the particle experiences a negative axial 
displacement (cf. Fig. flBTb)). Previous authors ([27]) also reported that even for Reynolds 
number 10, at certain points near the axis , the longitudinal displacement of the fluid particle 
can be negative. Fig. HBT c) shows that the particle period is nearly equal to the peristaltic wave 
period for a shear thinning fluid with n = 2/3. A comparison of Fig. IToT d) with Fig. fl5T c) and 
Fig. fTBT b) reveals that within one wave period, in the case of a shear thickening fluid, a fluid 
particle traverses more distance than in the case of a shear thinning/Newtonian fluid. 

4 Validation of the Results 

The last paragaph of Sec. 2 serves a validation of the analytical expressions derived in the present 
analysis. In order to validate the numerical results of the present study, we have compared our 
numerical results with those reported earlier by Takabatake and Ayukawa [27J on the basis of 
their numerical study of two-dimensional peristaltic flows. The comparison given in Fig. [7(a) 
shows that our results are in good agreement with those given in [27]. More particularly, one 
may observe from Table 1 that our results presented in Figs. [2jb,c) for the axial velocity of 
blood in peristaltic movement match well with the corresponding results reported in [27]. Also, 
a comparison between Figs. H](b) and HJ^e) shows that our results for the velocity contour are 
similar to those presented earlier by Selvarajan et al. [57]. Moreover, as shown in Figs. [2(a) the 
axial velocity variation in our case is in perfect matching with that reported by Shapiro et al. 
[TB] . Lastly, our results presented in Fig. ITBTa) are in excellent agreement with those reported 
in PS]. 

5 Summary and Conclusion 

In this study, an attempt has been made to investigate the peristaltic motion of blood in the 
micro-circulatory system, by taking into account the non-uniform geometry of the arterioles and 
venules. Treating blood as a Herschel-Bulkley fluid, the effect of amplitude ratio, mean pressure 
gradient, yield stress and the rheological fluid index n on distribution of the velocity and wall 
shear stress, pumping phenomena, streamline patterns and pathlines are examined under the 
purview of the lubrication theory. The salient observations are as follows: 
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Figure 15: Trajectories of massless particles for Newtonian/non-Newtonian fluids at different 

locations (Q = 0, = 0.5, r = 0, A = 0); for Q=0; • • • for Q = 30 2 /(2 + 2 ); • initial 
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(i) At any instant of time, there is a retrograde flow region for both Newtonian and non- 
Newtonian fluids, when AP = and also for some negative values of AP. 

(ii) With an increase in the value of n and 0, the regions of forward/retrograde flow advance 
at a faster rate. 

(iii) The parabolic nature of the velocity profiles are significantly affected by the value of 
the rheological fluid index 'n'. 

(iv) In case of a shear thinning fluid, flow reversal does not occur in the case of a channel 
having uniform geometry; it transforms to forward flow in the case of a diverging channel. For a 
shear thickening fluid, flow reversal reduces, but it does not vanish altogether when AP changes 
from to —1. 

(v) Non-uniform geometry affects quite significantly the distribution of velocity and wall 
shear stress as well as the pumping phenomena and other flow characteristics. 

(vi) Peristaltic pumping characteristics as well as the distribution of velocity and wall shear 
tress are strongly influenced by the amplitude ratio <fi and the rheological fluid index 'n'. 

The study bears the potential of significant application to the field of biomedical engineering 
and technology, because it is known that in roller pumps, the fluid elements are quite prone to 
significant damage of fluid elements; also, in the process of transportation of fluids in living 
organisms, by using arthro-pumps, the fluid particles are likely to be appreciably damaged. 
Moreover, the qualitative and quantitative aspects of the present investigation for the wall shear 
stress have a significant bearing on extra-corporeal circulation, where the heart-lung machine is 
usually used. In this case, there is a possibility that the erythrocytes of blood may get damaged. 
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